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Employing massive antennas array at the user terminals can be feasible by 
using millimeter wave (nmmWave) transmission which significantly reduce 
the antennas array size. The implementation of massive multiple input 
multiple output (MIMO) at the user terminals facilitates accurate 
beamforming. In this paper, a modified orthogonal matching pursuit (OMP) 
algorithm is used to design a wideband hybrid combiner based on the sparse 
structure of mmWave channel and orthogonal frequency-division 
multiplexing (OFDM). Based on OFDM, the wideband channel considered 
as multiple narrowband channels so a modified narrowband hybrid combiner 
can be implemented for each subcarrier channel in a manner where the RF 
combiner is the same for all subcarriers, whilst the baseband combiner is 
obtained for each subcarrier. For the multiuser 5G system, a wideband 
hybrid precoder based on the block diagonalization (BD) method is used at 
the base station (BS) to cancel the interference at each user due to the other 
users. The performance of this hybrid beamformers (precoder/combiner) are 


tested for different scenarios of base station antennas number, numbers of 
users’ antennas, and number of users. 
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1. INTRODUCTION 

The main technology of the fifth generation (5G) system is the large-scale antenna array. Massive 
multiple input multiple output (MIMO) system is a large-scale antennas array with huge antennas numbers 
up to several hundreds or thousands [1], [2]. Massive MIMO system increases the system capacity and 
provide more reliability to communication link leading to enhance energy and spectral efficiencies of the 
system. In massive MIMO, a huge number of antennas can be implemented at base station (BS) [3]-[5]. 
However, it is difficult to apply large number of antennas at the user terminals due to the large physical size 
of antenna. For example, the size of (4x4) uniform planar array (UPA) operates at a 5 GHz frequency is 
about 80 cm?. Whilst using millimeter-wave transmission (mmWave) leads to more realistic implementation 
of the massive MIMO technology at the user terminals [6], [7]. It can be seen that the size of the UPA 
mentioned above cab be reduced to about 2.6 cm? when the operating frequency is 28 GHz, even so about 
(16x16) UPA can be implemented in the 80 cm? area when operates at s28 GHz which provide more accurate 
beamforming [8], [9]. 

Analog beamformer is a beamformer implemented at the radio frequency (RF) stage, while digital 
beamformer implemented at the baseband stage [10]. The analog beamformer uses phase shifters to aim the 
main beam toward the desired direction [10], [11]. The beamforming process in the digital beamformer is 
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accomplished by means of the digital signal processing. It is kindly easy to do the beamforming process 
using digital beamformer with more accuracy and flexibility rather than analog beamformer [12]. However, a 
dedicated RF chain is required for each antenna element in digital beamforming. This RF chain module 
contains analog-to-digital converter (ADC), digital-to-analog converter (DAC), power amplifier, and mixers 
which makes its design at mmWave complex and expensive and large number of antennas are 
used [13], [14]. 

A structure of both digital and analog beamformers was introduced to take advantage of the both 
types of beamforming, and this combination called hybrid beamforming [15], [16]. Hybrid beamforming was 
mentioned for the first time in [17]. In hybrid beamforming, the antennas of array divided into many 
subarrays each connected to one RF chain only. So, the number of RF chains is significantly reduced, leading 
to a cheap and a simple system design. The performance gap between fully digital and hybrid beamformers 
can be reduced by optimizely design the hybrid beamformer. Also if the number of chains is more than two 
times the number of data, the gap becomes negligible [10], [18]. 

The amplitude and phase of the transmitted symbols can be modified using the digital beamformer, 
but the analog beamformer can modify only the phase of the transmitted symbols because it has a constant 
amplitude 1 IV (N+). In fact, the analog beamformer has quantized phases of N bits resolution, so accurate 
analog beamforming required a high resolution phase shifter which needed more bits [19], [20]. In this paper, 
a new design of wideband hybrid beamformer based on orthogonal frequency-division multiplexing (OFDM) 
is proposed for a massive MIMO multiuser 5G system. Since the mmWave channel is specular and has a 
sparse structure, the beamforming design problem can be formulate as a sparse reconstruction problem [21]. 
a modified algorithm of orthogonal matching pursuit (OMP) algorithm is used to find the optimal wideband 
hybrid combiner. The wideband 5G channel is divided into multiple narrowband channels by means of 
OFDM technique, so the proposed combiner can be applied to each subchannel. To eliminate the 
interferences between users, block diagonalization (BD) is used at the BS as a wideband hybrid precoder 
according to the modified joint spatial division multiplexing (JSDM) design mentioned in [12]. The design 
separation between precoder and combiner is accomplished giving Kronecker channel model [22], [23]. 

The following notation are used here: (.)~1,(.) and (.)¥ denote the inverse, conjugate and 
conjugate transpose respectively; || ||zand || || are -norm and Frobenius norm respectively. Statistical 
expectation is denoted by E[.]; diag(.) is the diagonal elements; CN (m, M) is a complex Gaussian 
distribution with mean m and covariance matrix M; span(V) and span+ (V) are the column space of (V) 
and its orthogonal complement respectively. SVD is the singular value decomposition and EVD is the eigen 
value decomposition. 


2. SYSTEM MODEL 
Consider a downlink multiuser MIMO-OFDM 5G system, shown in Figure 1, operates at mmWave 
frequency 28 GHz, where BS supplied with N, antennas placed in an UPA. The BS communicates with N, 


users through N, independent data streams such that N, = aH N; i, and N, ; represents the streams per user. 


Each user equipped with N,, antennas so N, = we Nya. Nar Neva» Nae and Ng, are the total number of 
RF chain at the BS, the number of dedicated RF chain at the BS for the i‘” user, the total number of the 
receiver RF chain of all users, and the number of receiver RF chain for it” user respectively. Working at the 
mmWave band can enhance the beamforming and transmission bit rate, but to obtain these benefits, OFDM 
is needed to deal with the large wideband of mmWave channel by dividing this wideband channel into 
several narrowband channels, so a conventional narrowband beamformer can be applied to this 
system [12], [24]. 

The hybrid precoder is applied at the BS, while the hybrid combiner is applied at the user terminals. 
At both side of the transmission link, the analog precoder/combiner is common over all subcarrier channels, 
and the digital beamformer is applied to each subchannel. f¥ and w* are the hybrid subcarrier precoder and 
hybrid subcarrier combiner of the kt” subcarrier channel for it user respectively. The f* of the kt” 
subchannel is of the form Fr, Tees where analog precoder Fre; E CNEXNRF A, and baseband precoder 
fasa E CXRFa*Nsi. The w¥ is of the form WrriWbe i> where Wrz i E CNr*NRrF4 ig the analog combiner, and 
Whi E CNkri*Nsi is the digital combiner. 

In this paper, the focus is on the design of the hybrid subcarrier combiner while adopting the design 
mentioned in [12] for the hybrid subcarrier precoder. In fact, the design of the hybrid subcarrier 
precoder/combiner can be separable as will be showed later. The digital subcarrier precoder fay maps the 
transmitted symbol signal x¥ of kt subcarrier and belongs to 4t" user to the input of OFDM transmitter, 
then Nerd RF chains is used to upconvert the OFDM signal to RF stage. Finally, the output of the chains is 
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connected to the BS antennas array using Néri x N; phase shifters of the analog precoder Fr, because of 
using the fully connected structure. 
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Figure 1. Multi-stream multi-user massive MIMO-OFDM 5G system 


So, the transmitted signal y* as (1): 
yi = FrrifpeiXt (1) 


The received signal, before the hybrid subcarrier combiner, rj‘ at i” user as (2): 


= JP; HY Fer afeeixt + Dp, VP; Al Fer. fio xf +n (2) 


jti 


The noise signal n; has complex gaussian distribution e.g., n; ~ CN (0, Rn), Rn; is the 
covariance matrix, and the average power at the 4t? user is P;. The wideband channel Hf e CNri*Ne 
between BS and the it? user at k OFDM subcarrier [25]. Following the approach of [12], the MIMO 
channel HE becomes as (3): 


Nei 


= VaLa” Lm Bess rd Pimi Ohmi) Dima Oimai) (3) 
This channel model consists of Nop; clusters, and Np; rays per t? cluster, so Nays, ay Nei- 
Yi = JNLN,/Nraysi> and Bema is the complex gain of the m*” ray in the #*” cluster at the k*” subchannel, 
whereas Pima ema) is the azimuth (elevation) angle of arrival whilst Ce is the azimuth 
(elevation) angle of departure. a,(.) and a,,(.) are the antennas array response vectors of BS and the i” 
user respectively. The right hand of (3) is similar to the narrowband channel model [26], so a narrowband 
beamformer can applied for each OFDM subcarrier. For simplicity, the following assumption are made: each 
user has the same number of data streams, so {P; = P,W i}, x¥ has a normal distribution e.g., 


: x 2 A f 
~ N (0, Rxi), the covariance matrix Rx; = 1/y ess { [Far ifto.al, = NY k, ij for the transmitted 
sa Ns ; / 
power constraint. Assume r¥ propagates through a block fading narrowband channel, this assumption is true 


as long as using OFDM technique [24], [27]. Moreover, based on the Kronecker channel model [22], [28] the 
hybrid subcarrier precoder and combiner can design separately. 


3. HYBRID BEAMFORMING DESIGN 
The essential key in the design of hybrid beamforming is finding the optimal precoder and combiner 
weights fisi FRF wkp Wrr, that minimize the mean square error (MSE) between transmitted and 
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received symbols for all users and over all subcarriers. As mentioned earlier it is acceptable to firstly design 
the hybrid subcarrier precoder then the hybrid subcarrier combiner. The design of hybrid precoder and 
combiner is described in the following sub-sections. 


3.1. Hybrid subcarrier precoder design 

The JSDM algorithm [29] and its modification [12] show a little performance degradation in 
comparison to the optimal fully digital precoder and reduce the overhead information needed to fed back, so 
this algorithm is adopted here to design the hybrid subcarrier precoder for the system in Figure 1. The users 
assumed to be distributed in a manner that they are well-spatially separated and there is only one user per 
group. The powerful of this algorithm is that the analog precoder design depends on the second-order 
statistics of the channel and not on its first order statistics. The analog precoder Fpp; is designed according to 
the procedures listed in Algorithm 1 [12]. 


Algorithm 1: design of analog precoder Fprri 
1. Estimate HŽ. 
2. Define M = E[H¥ 4HF]. 
3. Find the eigenvectors Vë by computing EVD(M*). 
4 Feedback V# to the BS. 
5. Computing V; by averaging vk over all subchannels. 
6 
7 
8 


Define £,=[V,V, Va Vier ~ Vnu Vm) - 

Find the eigenvectors B; by computing SVD(E;). 
Define B, such that Span (8) = Spant ({E;: vj # i}). 
9. BS broadcasts BO for all users. 

10. Define HK = HG, 

11. Define ME = E[HE "HE ]. 

12. Find the eigenvectors Vf by computing EVD(M;) 

13. Computing y, by averaging vE over all subchannels. 
14. Find ÑP such that ¥=[9,2 y,®]. 


yxa 
15. Finally, the analog precoder is defined, Foz, = 8,” * yo 


: : x : i oe H 

Now, the channel is defined as H* = HE Fras and its covariance matrix is Wk = Frei ME Fara. 

Since, Feri" Fer, = Cl Niri for large number of BS antennas [21], [24], the digital subcarrier precoder fp; 
; ; 4 ; 


is the normalized N,; eigenvectors of M, HA 


3.2. Hybrid subcarrier combiner design 

Since the hybrid subcarrier combiner design is independent for each user but it is the same for all 
users. Therefore, the procedures of the design discussed here are for it” user only. The second term of the 
right-hand side of (2) is eliminated because the column vectors of Fpp; eliminates inter-group interference 
between users [30], and (2) as (4): 


rë = JP, H} Freafopaxe + n (4) 


By applying the hybrid subcarrier combiner w*, the received signal 7# in (5): 


k_ kik 
_ k H H prk k „k k H H 
= 4P, wri Wprri HiFrrifBBiXi + Bei Wrri ni (5) 


The hybrid subcarrier combiner is designed to minimize the MSE between xf, and rj. The 
optimization problem is written as (6): 


(Wer ogee) = argmin Elljxk - rE] 
WRF i» OBB i 


S.t. Wpgra © Wer (6) 


t _ kopt ; ‘ mae ‘ i 
Where Wot o Opay are the optimal values of the analog combiner and digital subcarrier combiner 


respectively. Wp, is the feasible set of analog combiners. The optimal unconstrained solution wri P of (6) is 
well known [31] as (7): 
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wE” = [efe tre lB epe] 


-174 
VPI H[ -Pi H 
=|= föga Frea” Hi (Fit Fans fofb. HFera He + shat. 
Ng ji Ng i É 
eF H 
of AN si H 
ole =( fis. "Frea" HE Hi Porafloa + ly.) flea Frea" HE | (7) 


Where (a) is obtained by applying the push-through identity of matrix inversion lemmas. In order to apply 
this optimal combiner dal to the system of Figure 1, this combiner needs to be decomposed into the 
product of analog and digital subcarrier combiner Wor. ORE pi with Wer; E Were. Due to this complex 


non-convex constraint, there is no optimal solution of (6). Since each subcarrier channel is modelled as a 
narrowband channel, the approach of [21] is applied here to overcome this difficulty. According to this 
approach the optimization problem can be rewritten as (8). 


1 
opt | k,opt\ _ z k „kH |2 k,Opt k 
(Wri ® Opp i ) = argmin ||E |r ri (w; — Wrri wks) 


k 
WRF i» WBB i F 


S.t. Wpgra © Wer (8) 


Now the optimization problem is to find the E [rE rk” |-weigntea Frobenius norm projection of the 


unconstrained combiner wirt onto the set of hybrid subcarrier combiners of the form Wp- i wkp with 
Wrri E Ware. Again, the feasibility constrain makes the solving of (8) very difficult. Fortunately, due to the 
sparse structure of H*, a near optimal combiner can be found by more restricting Wer; So its columns are 
Neri columns of an N,i*Nraysi matrix of receiver array response vectors 


Ari = EEU Oraa) ari Cer ee || and solving as (9). 


1 
~k,opt _ y k kH |2 k,Opt <k 
BBi = arg min E |r ri | (w; — Ár i Dep a) 
OBBi F 
s.t. [diag ( whe w wep “I, = Neri (9) 


Algorithm 2: design of hybrid subcarrier combiner wk 
A. for k < number of OFDM subcarrier do 
a. Initialization 
1. Estimate H* 
2. Find fg, and Fer; according to Algorithm 1 
3. Compute wE according to (7) 
4. Let Wiri be an empty matrix 
5. Compute Er = E [rr] 
6. Set Wres = w” 
b. Construction of Whe i and Wiri 
1. for J = 1,2, Ngri do 
2. D=A GEL Wrest? 
3. £=argmax(diag(bo*)) 
4. Add £ column of Api to Wri 
5 


; ; H = H k 
Using least squares (LS) to find wk = (Wira ek whee) Wir Eka r 
6 w D= ww PW OBB A 
* = k, 
a Je; Se -Wiriokeal p 
7. end for 
B.end for 


C. Construction of Wprp; by taking the mean of Wiri over all subcarriers 


Now the problem is equivalent to sparse signal recovery problem and can be solved through OMP 
[21]. Ok, i € CNraysi*Ns acts as an auxiliary variable from which the digital subcarrier combiner is obtained. 


mk — Nr =k r 
The constrain ||diag (a bei Ope i DIN = Neri ensure that the @;, ; cannot have more than Npr; non-zero 
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rows which form the digital subcarrier combiner Win ds thus only Ng, columns of A, is effectively selected 


to construct the analog combiner Wpr ¿. The complete design of hybrid subcarrier combiner w* is illustrated 
in Algorithm 2. 


4. RESULTS AND DISCUSSION 

In this section, the performance of Algorithms 1 and 2 of the hybrid subcarrier beamformers 
fE and w design is evaluated considering the system in Figure 1. Different numbers of users randomly 
distributed within a one-thousand-meter diameter are used to test the performance of these algorithms. To 
ease the calculations, all users are equipped with the same number of antennas and communicates two data 
streams Ns; = 2 with BS. The wideband 28 GHz channel is dividing into 256 OFDM subchannels plus 64 
cyclic-prefix (CP). The clustered channel in (3) is adopted to model these subchannels, so each subchannel 
consist of ten clusters each with ten rays. For channel estimation a channel sounding system is used. The 
results of proposed algorithm are compared to that of fully digital beamformer for different numbers of BS 
array antennas, also different numbers of users’ array antennas are used to test the system in Figure 1. The 
antennas of both BS and users’ array are arranged in UPA, and the number of transmitted RF chain is 3 per 
data stream, Neri = 3 x Ni and the receiver RF chain is one per data stream, Nap; = Ns. 

First, a comparison between a fully digital beamformer with full channel state information 
(FCSIFD) and the proposed beamformer which designed based on Algorithms 1 and 2 is presented in 
Figure 2 in term of average spectral efficiency versus signal-to-noise ratio (SNR). The comparison is done 
for a various number of BS antennas N,=64, 128 and 256, and the number of users is six. In Figure 2(a) The 
number of antennas per user is four, whilst is six for Figure 2(b). The performance of the proposed 
beamformer (Algorithms 1 and 2) is also compared with that of the precoder designed using Algorithm 1 to 
show the enhancement obtained by applying the combiner at the users’ terminals. 

As shown in Figure 2, increasing the number of BS antennas reduces the gap between FCSIFD and 
Algorithms 1 and 2. Also, it has been shown that the performance of Algorithms 1 and 2 superiors to that of 
Algorithm 1 especially when the number of users’ antennas increases which means that the applying of 
combiner at users will enhance the system in Figure 1. So, the gap between Algorithm 1, Algorithms 1 and 2 
is about 0.97 at SNR=10 dB for N, = 64 antennas and N,; = 4 antenna, but the gap increases to be 
approximately 1.5 when the number of users’ antennas is increased to be 8 antennas. Moreover, the 
performance gap between Algorithms 1, 2 and Algorithm 1 is more when the number of BS antennas is low, 
e.g., Ne = 64, so it is found that the gap is about 1.5 at SNR=10 dB for N, = 64 and N,.; = 8 whilst the gap 
is about 1.2 when N, = 256. It means that the proposed Algorithms 1 and 2 provides more spectral 
efficiency for system with low number of BS antennas. 
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Figure 2. Spectral efficiency of FCSIFD, Algorithm 1, Algorithms 1 and 2 for six users, each user equipped 
with; (a) four antennas, and (b) eight antennas, two data stream per user 


Figure 3 shows the performance of FCSIFD, Algorithm 1, Algorithms 1 and 2 for the system when 
N, = 6 users, and N,; = 16 antennas for 128, and 256 BS antennas. The system with 64 BS antennas can’t 
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be implemented here because the total number of users’ antennas, N, = 96, is more than that of BS. Here, 
the performance gap is about 1.89 when N, = 128, while it is about 1.2 when N, = 256. Again, it has been 
shown that increasing the number of users’ antennas enhances the performance of Algorithms 1 and 2 and 
reduce the gap with FCSIFD. To study the effect of the number of users’ antennas on the system in Figure 1, 
the spectral efficiency versus number of antennas per user, N, į = 4,8, 16, and 32, is plotted in Figure 4 for 
different number of users and BS antennas where the SNR=5 dB. 

The results of Figure 4 prove that the performance of Algorithms | and 2 is superior to that of 
Algorithm 1 specifically when increasing the number of antennas per user. The enhancement of Algorithms 1 
and 2 is about 1.16 above Algorithm 1 when N; = 256, N, = 6, and N,; = 8, Correspondingly the 
enhancement increases to about 1.4 when N,; = 16, and exceeds 2 for N,.; = 32. Again, the system with 
N, = 128 cannot implemented when the number of antennas per user N,.; = 32, because the total number of 
users’ antennas N, = 192 which is much more than the number of BS antennas. 
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Figure 3. Spectral efficiency of FCSIFD, Algorithm 1, Algorithms 1 and 2 for six users; each user equipped 
with sixteen antennas, two data stream per user 
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Figure 4. Spectral efficiency of Algorithm 1, Algorithms 1 and 2 versus number of antennas per user for six 
and eight users 


Moreover, it shown that the performance of Algorithms 1 and 2 when N,, = 8 is approximately the 
same as that of Algorithm 1 when N,, = 6 for either 128 or 256 BS antennas despite that the performance 
reduces as the number of user increases. The relation between spectral efficiency and the number of users is 
illustrated in Figure 5. Here, it can be shown that the performance of the system in Figure 1 decreases as the 
number of users increases. However, the performance of Algorithms 1 and 2 when N,.; = 8 is almost the 
same as the performance of Algorithm 1 when N, ; = 16 for both 128 and 256 BS antennas. This means that 
applying the proposed combiner designed based on Algorithm 2 either enhances the system or decreases the 
number of users’ antennas required to achieve a certain sum rate. 
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Figure 5. Spectral efficiency of Algorithm 1, Algorithm 1 and 2 versus number of users 


5. CONCLUSION 

In this paper, an OFDM-based wideband hybrid subcarrier beamformer is designed for multi users 
MIMO 5G system. The design consists of two stages. In the first stage, the hybrid subcarrier precoder is 
designed based on the concepts of BD, whereas the hybrid subcarrier combiner is designed based on the 
concept of sparse signal recovery. Both the precoder and combiner is designed for wideband mmWave 
system through using OFDM technique to deal with this wideband channel as multiple narrowband channels. 
This wideband beamformer has been tested for a massive MIMO 5G system where the effects of the number 
of BS antennas, the number of antennas per user and the number of users on the performance of the proposed 
beamformer are evaluated. The results show that the proposed beamformer has an acceptable degradation in 
spectral efficiency in comparison to the performance of the fully digital beamformer with full CSI. Also, they 
show that the proposed beamformer (combiner plus precoder) provides good enhancement in comparison to a 
system applying precoder only where the enhancement of the spectral efficiency exceeds two when the 
number of antennas per user is more than 32. Despite the proposed beamformer is not an optimal beamformer 
but it provides a good performance in term of spectral efficiency and computational efficiency because it 
depends on the second channel statistics which are change slowly with time, so this eventually reduces the 
overhead information need to deal with them. 
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